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Abstract 

We point out that determination of the MNS matrix element \Ue 3 \ = S 13 
in long-baseline ^ Ve neutrino oscillation experiments suffers from large 
intrinsic uncertainty due to the unknown CP violating phase <5 and sign of 
Am^g. We propose a new strategy for accurate determination of 9iy, tune the 
beam energy at the oscillation maximum and do the measurement both in 
neutrino and antineutrino channels. We show that it automatically resolves 
the problem of parameter ambiguities which involves <5, 613 , and the sign of 
Amfg. 
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I. INTRODUCTION 


With the accumulating evidences for neutrino oscillation in the atmospheric [Q, the 
solar [§] and the accelerator neutrino experiments it is now one of the most important 
subjects in particle physics to explore the full structure of neutrino masses and the lepton 
flavor mixing. In particular, it is the challenging task to explore the relatively unknown (1-3) 
sector of the MNS matrix namely, ^ 13 , the sign of Am ^3 and the CP violating phase 
6 . The only available informations to date are the upper bound on ^13 from the reactor 
experiments |^, and an indication for positive sign of Am^g by neutrinos from supernova 
1987A [p. Throughout this paper, we use the standard notation of the three flavor MNS 
matrix, in particular = 8136 “*^, and define the neutrino mass-squared difference as 


Am'fj = rrij — rrif. 

The long baseline z/g neutrino oscillation experiment is one of the most promising 

way of measuring ^ig. In particular, it is expected that the JHF-Kamioka project which 
utilizes low energy superbeam can go down to the sensitivity sin^ 2^13 ~ 6 x 10“^ |^. A 
similar sensitivity is expected for the proposed CERN —Frejus experiment [§]. Although a 
far better sensitivity is expected to be achieved in neutrino factories [^, it is likely that the 
low energy conventional superbeam experiments are the ones which can start much earlier. 
Therefore, it is of great importance to examine how accurately ^13 can be determined in this 
type of experiments. 

In this paper, we point out that determination of sin^ 2^13 by using only neutrino channel 
suffers from large intrinsic uncertainty of ± (30-70) % level due to the unknown CP violating 
phase 5 and the undetermined sign of Amfg. It should be noted that the intrinsic uncertainty 
exists on top of the usual experimental (statistical and systematic) errors. To overcome the 
problem of the intrinsic uncertainty, we suggest a new strategy for determination of ^13 by 
doing appearance experiments utilizing both antineutrino and neutrino beams. Our proposal 
is a very simple one at least at the conceptual level; tune the beam energy to the oscillation 
maximum and run the appearance experiments in both —>■ z/g and z/^ —>■ Ve channels. 

We will show that it not only solves the problem of intrinsic uncertainty mentioned 


above but also resolves the (5— 6 * 13 ) two-fold ambiguity discussed in Ref. [|T3- Furthermore, 
it does not suffer from possible ambiguity due to the unknown sign of Am^g, the problem 
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first addressed in Refs. Hu-Ill We are aware that there are combined ambiguities to be 
resolved (even ignoring experimental uncertainties) to determine a complete set of param¬ 
eters including S, di 3 , and the sign of Am^g, which are as large as four-fold We take 
experimentalists’ approach to the ambiguity problem and try to resolve them one by one, 
rather than developing mathematical framework for the simultaneous solutions. The most 
important issue here is again to accurately determine ^ 13 , because then all the combined 
ambiguities will be automatically resolved, as we will show below. 

II. INTRINSIC UNCERTAINTY IN DETERMINATION OF DUE TO CP 

VIOLATING PHASE 

Let us clarify how large uncertainty is expected for determination of 6^13 due to our 
ignorance of 6 in the > z^e appearance experiment. To achieve intuitive understanding of 
the issue we use the CP trajectory diagram introduced in previous papers ||TT|JT^ . Plotted in 
Fig. 1 are the CP trajectory diagrams in bi-probability space spanned by P^i^) = t'e) 

and T’(h) = ^ he) averaged over Gaussian distribution (see next paragraph) with three 

values of ^ 13 , sin^ 26*13 = 0.05 and 0.02 for > 0 case and sin^ 26*13 = 0.064 for Am^g < 0 
case. Since we assume lAm^gl 3> |Am 32 | the sign of Am^g is identical with that of Arngg. 
(The fourth one with sin^ 26*13 = 0.04 is for our later use.) The values of sin^ 2^13 for the 
second and the third trajectories are chosen so that the maximum (minimum) value of (P(z/)) 
of the second (third) trajectory coincides with about 1.1 %, the minimum value of {P{u)) 
of the hrst trajectory. The remaining mixing parameters are taken as the best £t value of 
the Super-Kamiokande (SK) and the K2K experiments [|^, lAm^gl = Am^^^ = 3 x 10“^ 
eV^, and the typical ones for the large mixing angle (LMA) MSW solar neutrino solution as 
given in the caption of Fig. 1. 

While we focus in this paper on the JHF experiment with baseline length of 295 km, 
JAERI-Kamioka distance, many of the qualitative features of our results remains valid 
also for the CERN-Frejus experiment. Throughout this paper we take the neutrino energy 


Our new strategy and these results were announced in the 8th Tokutei-RCCN workshop 
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distribution of Gaussian form with width of 20 % of the peak energy. Of course, it does not 
represent in any quantitatively accurate manner the effects of realistic beam energy spread 
and the energy dependent cross sections. But we feel that it is sufficient to make the point 
of this paper clear, illuminating our new strategy toward accurate determination of ^ 13 . 

Suppose that a measurement of appearance events gives us the value of the oscillation 
probability {P{v)) ~ 1.1 %. Then, it is obvious from Fig. 1 that a full range of values of 
sin^26'i3 from 0.02 to 0.064 are allowed (even if we ignore experimental errors) due to our 
ignorance to the CP violating phase 5 and the sign of If we know that the sign is 

positive, for example, the uncertainty region would be limited to 0.02-0.05, which is still 
large. 

Let us estimate in a syatematic way the uncertainty in the determination of ^13 due to 
the CP violating phase 5. To do this we rely on perturbative formulae of the oscillation 
probabilities P’(z^) and P{V) which are valid to first order in the matter effect |^. With 
relatively short baseline ~ 300 km or less the first-order formula gives reasonably accurate 
results. Ignoring 0(sin^26'i3) terms the formula can be written with use of the notation 


Ajj — 


Am?. L 


- ».7 


2E 


{L and E denote baseline length and neutrino energy, respectively) in the form 


Pi^v/i/) = sin^ 26^13 + 2 Q sin 26^13 cos ^ 


A 


13 


±5 


( 1 ) 
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P±(Ai 3) = s 


sm 


A 
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( 2 ) 


Q = C12S12C23S23 Ai 2 sin 


A 


13 


( 3 ) 


^ It may be worth to remark the following: Low energy neutrino oscillation experiments with 
superbeams are primarily motivated as a result of the search for the place where CP violating 


effects are comparatively large and easiest to measure [|^. See e.g., |T^ for works preceding to 


|]T|. Unfortunately, this large effect of 6 is the very origin of the above mentioned large intrinsic 
uncertainty in determination of 613 . 
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where a = \/2GpNf. denotes the index of refraction in matter with Gp being the Fermi 
constant and N^. a constant electron number density in the earth. The ± signs in P± refer 
to the neutrino and the antineutrino channels, respectively. 

The maximum and the minimum of P{i') for given mixing paramters, neutrino energy 
and baseline is obtained at cos (^5 + = ±1. Then, the allowed region of sin 2^13 for a 

given value of P{i^), assuming blindness to the sign of Am^g, is given by 


- \Q\ ^ ^ + p+(-Ai3)pm + 101 

-- \ Sin \ -- 


(4) 


P+(Ai3) - P+(-Ai3J 

In Fig. 2 presented is the allowed region of sin^ 2^13 for a given value of measured oscillation 
probability P{i^). Figures (a)-(c) correspond respectively to neutrino energies (a) E = 500 
MeV, (b) 716 MeV (oscillation maximum), and (c) 1 GeV. One notices that the intrinsic 
uncertainty is large. It strongly depends on the value of sin^ 2^13 and gradually decreases 
as E grows. Roughly speaking, it ranges between, 45 % (30 %) at sin^ 26*13 = 0.1 and 
~ 80 % (70 %) at sin^ 26*13 = 0.01 at = 500 MeV (716 MeV). Notice that all the results 
shown in the plots in this paper were obtained by numerically solving the neutrino evolution 
equation assuming constant matter density without using the hrst-order formula. 

The size of the intrinsic uncertainty must be compared with the statistical and the 
systematic errors which are expected in the actual experiments. A detailed estimation of the 
experimental uncertainties is performed for the JHF experiment by Obayashi assuming 
the off-axis beam (OA 2 ) and running of 5 years. The results strongly depend upon 6 * 13 . 
We quote the case of three typical values; sin^26*i3 = 0.1-1-0.018 — 0.014, 0.03-1-0.010 — 0.007, 
and 0.01 -|- 0.007 — 0.006. The errors include not only statistical but also systematic ones. 
We implemented these errors in Fig. 2b which is drawn with the similar energy as the 
peak energy of OA2 beam (~ 780 MeV). We should note, however, an important difference 
between Fig. 2 and the plot in [^; the abscissa of Fig. 2 is the Gaussian averaged probability, 
whereas the corresponding axis of the plot in is the number of events. Therefore, we 
tentatively determined the location of errors in Fig. 2 so that the center of the error bars 
coincide with the center of the allowed band of sin^ 26 * 13 . Keeping this difference in mind, 
we still feel it informative for the readers to display the expected experimental uncertainties 
in Fig. 2b for comparison. 
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Therefore, the intrinsic nncertainty dne to S and nndetermined sign of is larger than 
the expected experimental errors in most of the sensitivity region for ^13 in the experiment. 
We note that the experimental errors are dominated by the statistical one in phase I of 
the JHF-SK nentrino project and hence it should be improved by a factor of ~ 10 in two 
years of running in the phase II with a megaton water Cherenkov detector |Q. Thus, the 
intrinsic uncertainty completely dominates over the experimental ones if one stays only on 
the neutrino channel. 

III. POSSIBLE WAY OUT AND THE RELATIONSHIP WITH 6»i3 - J 

AMBIGUITY 

Let us discuss possible ways out of the uncertainty problem in the determination of ^ 13 . 
It is tempting to think about seeking better resolution by adding more informations. A 
natural candidate for such possibilities in this line of thought is to do additional appearance 
experiment > Pe using antineutrino beam. While it strengthens constraints, it does not 
completely solve the uncertainty problem even if we ignore the experimental errors. It is 
due to the inherent two-fold ambiguity which exists in simultaneous determination of 6 and 
6*13 us has been pointed out by Burguet-Castell et ah 0 . While their discussion anticipates 
applications to neutrino factory, the issue of the two-fold ambiguity is in fact even more 
relevant to our case because of the large effect of 6 as we saw in the previous section. 

The existence of the two-fold ( 6*13 — 5) ambiguity is easy to recognize by using the CP 
trajectory diagram. We show in Fig. 1 by a dash-dotted curve another trajectory drawn 
with sin^ 20 i 3 = 0.04 which has two intersection points with the solid curve trajectory 
with sin^ 26*13 = 0.05. Suppose that measurements of neutrino and antineutrino oscillation 
probabilities P{i') and P{v) had resulted into either one of the two intersection points. 
Then, it is clear that we have two solutions, for positive Am^g, (sin^ 26 * 13 , <^) = (0.04, O.OStt) 
and (0.05 , 0 . 357 r) for the upper intersection point, and (sin^ 26 * 13 , = (0.04, 1 . 47 r) and (0.05, 

I.Ttt) for the lower intersection point. Similar two-fold ( 6*13 — S) ambiguity also exists for 
negative Ami 3 which however is not shown in Fig. 1. In other word, we can draw two 
different CP trajectories which pass through a point determined by given values of P{v) and 
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F(p). This is the simple pictorial understanding of the (^13 — two-fold ambiguity which 
is uncovered and analyzed in detail in |]^. We will show in the next two sections that the 
ambiguity is automatically resolved by our proposal. 


IV. NEW STRATEGY FOR DETERMINATION OF 6»i3 


We now present our new strategy for determination of ^13 which avoids the problem 
of the large intrinsic uncertainty. It is intuitively obvious from the CP trajectory diagram 
displayed in Fig. 1 that if one can tune the experimental parameters so that its radial 
thickness (which measures the cos 5 term in Eq. (|lD) vanishes then the two-fold ambiguity 
is completely resolved. It occurs if we tune the beam energy at the oscillation maximum so 
that Ai 3 = tt as is clear from Eq. ©• 

We explain below in more detail how it occurs and then discuss by what kind of quan¬ 
tity 013 is determined. In the following discussion we assume that the mixing parameters 
lAm^gl ~ |Am 23 | = Ami 2 = Am^, 023 , and 0 i 2 are accurately determined by the 

time of the experiment. It is not so unrealistic assumption in view of the array of exper¬ 
iments ongoing (SK, SNO, K2K, KamLAND), on schedule (Borexino, MINOS, OPERA), 
or in planned (JHF). For example, the uncertainty in measurement of 623 is expected to be 
(5(sin^ 2023) ~ 0.01 in JHF phase I 0. 

We note that the oscillation probabilities (J^) can be written as 


F(z/) = A cos J -|- F sin J -|- C+ 

F(p) = AcosJ — FsinJ-|-C_ (5) 


where A = Q cos f^j, B = —Qsin(-^), and C± = F±sin^ 20 i 3 in the present approxi¬ 


\ 2 ) 


mation. It is easy to show from this expression that CP trajectory diagram is elliptic in the 
approximation that we are working 0 - (In fact, it is the case for all the known perturbative 
formulae.) Given (^) it is simple to observe that the CP trajectory is a straight line at the 
oscillation maximum, A = 0 ; the equation obeyed by the oscillation probabilities is given as 
F(z/) -|- P{v) = C+ + C-. Moreover, the hrst order matter effect cancels in C+ + C-, leaving 
the vacuum peace of P±. Therefore, the slope of the straight-line CP trajectory is the same 
as that in vacuum, and the matter effects affects only on the maximum and the minimum 
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points of the straight line in P(z/) and P{p) coordinates. Thus, once a set of values of P{y) 
and T’(h) is given by the experiments, one can determine C+ + C- as the segment of the 
“CP straight line” in the diagram, and hence sin^ 26*13 fo which is proportional. 

Thus, measurement of T’(z^) and P{V) at the oscillation maximum implies determination 
of 013 without suffering from any uncertainties due to unknown value of 5 and the sign of 
Am?3. 

In Fig. 3 we present the thinnest trajectories with the tuned value of the energy E = 760 
MeV for L = 295 km (JAERI-Kamioka distance) with sin^ 20 i 3 = 0.05 and 0.02, by taking 
the other mixing parameters given in the caption of Fig. 1. The energy would he E = 716 
MeV if we sit on the oscillation maximum. It arises because the contributions from higher 
and lower energy parts around the peak energy do not completely cancel because of the extra 
1/E factor in the cos 5 term for symmetric Gaussian beam width. Thus, we need slightly 
higher energy to have the thinnest trajectory. It should be noted, however, that the feature 
highly depends upon the specihc beam shape, and will also be affected by the fact that the 
cross section has an extra approximately linear E dependence. 

The slightly different slope of the straight-line trajectories of positive and negative 
indicates the higher order matter effect. This effect must be (and can be) taken into account 
when one try to determine 0 i 3 following the method proposed above. 


V. COMMENTS ON THE RELATIONSHIP WITH 6*13 - < 5 - SIGN OF AM? 
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AMBIGUITIES 


We now show that the ( 0 i 3 — S) ambiguity is automatically resolved by tuning neutrino 
energy at the oscillation maximum. It must be the case because two straight-line trajectories 
with the same slope do not have intersection points. For our purpose, it suffices to work 
with oscillation probability at a hxed monochromatic beam energy because averaging over a 
hnite width complicates the formalism and may obscure the essence of the problem. It can 
be shown |T^ that the difference between the true (^is) and the false ( 0 ( 3 ) solutions of 0 i 3 
for a given set of P{i') and P{V) is given under the small 0 i 3 approximation by 


0 iq — 013 — 


sin 5 — 2 ; cos 5 2Q 


1 + ^2 


P--P, 


sm 


A 


13 


( 6 ) 





where 


F- + F-. 


z = 


F_ - P. 


tan 


+ 


A 


13 


(7) 


Hence, the difference vanishes at the oscillation maximum, A 13 = tt, which means z —> cxo. 
Thus, no (013 — (5) ambiguity exists at the oscillation maximum as expected. 

It should be emphasised that our strategy of tuning beam energy at the oscillation max¬ 
imum is not affected by the ambiguity correlated with the sign of Amfg which is discussed 


in Ref. It is because the matter effect split the straight-line CP trajectories of positive 
and negative Aruig toward the direction of the line itself in first order of the matter effect. 
The possible correction comes from higher order matter effect which is small in the relatively 
short baseline of the JHF (as well as the CERN —Frejus) experiment, as shown in Fig. 3. 
The effect can be easily taken care of in the actual determination of 0 i 3 . 


VI. CONCLUDING REMARKS 

In this paper, we proposed a new strategy for accurate determination of 0 i 3 without 
suffering from the intrinsic ambiguity due to unknown value of S. That is, tune the beam 
energy at the thinnest CP trajectory and do the measurement both in neutrino and antineu¬ 
trino channels. We have shown that our new strategy completely resolves the ambiguities 
in the determination of 0 i 3 due to S and due to the sign of Amig within the experimental 
accuracy attainable in such experiments. 

One of the proposal which could be extracted from the strategy described in this paper 
is a possibility of having beam as early as possible. It would be the promising option 
for the case of relatively large sin^ 20i3, say, within a factor of 2-3 smaller than the CHOOZ 
bound. In this case, the —>■ z/g appearance events can be easily established in a few years 
of running of next generation neutrino oscillation experiments. Then, the uncertainties in 
determination of 0 i 3 would be greatly decreased by switching to beam rather than just 
running with the beam. 

What would be the implication of our strategy to the determination of 67 The tuning of 
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beam energy at thinnest trajectory in fact also provides a good way of measnring The 
ambiguity (5 —> vr — 5), however, is unresolved and it would necessitate supplementary mea¬ 


surement either by using “fattest” trajectory conhguration or by second detector with 
different baseline distance |^. We should emphasize that once 6*13 is measured accurately 
there is no more intrinsic ambiguities in determination of 6. We have explicitly shown that 
( 5 - 6 * 13 ) ambiguity is resolved. The only ambiguity which would survive (from the viewpoint 
of determination of 5) would be the accidental one that arises in a correlated way (5 — sign 
of Am^g), which is nothing but the remnant of (5 —>■ tt — 5) degeneracy in vacuum [^. But 
it is also resolved by either one of the two second measurements mentioned above. 


Nore added: 


While this paper was being written, we became aware of the paper by Barger et ah |21 


whose results partially overlaps with ours. However, most of the ambiguities discussed in 
the paper will be gone once 6*13 is determined accurately, as we noted above. 
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purpose of having maximal CP-odd (sin J) term at the energy |^,^. 
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FIGURES 



<P(v)> 

FIG. 1. CP trajectory diagrams showing the contours of (P(z^)) = ^e)) and 

{P{v)) = {P{v^ Ve)) as a function of 6. The Gaussian energy distribution of neutrino beam with 
(P) = 0.5 GeV with width o" = 0.1 GeV is assumed and the baseline length is taken as L = 295 km. 
The mixing parameters are fixed to be = ±3 x 10“^ eV^, sin^ 2023 = 1 - 0 , Am \2 = 6.2 x 10 ^ 

eV^, tan^ 6*12 = 0.35. We take the matter density as p = 2.8 g/cm^ and the electron fraction as R. 
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<P(v)> 

FIG. 2. Allowed region of sin^ 20i3 is shown as a shaded strip for given values of t'e)) 

(given in %) assuming Gaussian energy distribution of neutrino beam centered at {E) = (a) 0.5, 
(b) 0.716, and (c) 1.0 GeV with 20 % width a of {E) for L = 295 km. If the sign of is 

known, the allowed region is within the solid (Am 23 > 0) dashed {Arn^^ < 0) lines. The 

other mixing parameters and the matter density are taken as in Fig. 1. 
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<P(v)> 



The beam profile, the mixing parameters and the matter density are taken as in Fig. 1. 




